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Abstract 

We present exact calculations of the zero-temperature partition function for the g-state 
Potts antiferromagnet (equivalently, the chromatic polynomial) for two families of arbitrarily 
long strip graphs of the square lattice with periodic boundary conditions in the transverse 
direction and (i) periodic (ii) twisted periodic boundary conditions in the longitudinal direc- 
tion, so that the strip graphs are embedded on a (i) torus (ii) Klein bottle. In the limit of 
infinite length, we calculate the exponent of the entropy, W{q), show it to be the same for 
(i) and (ii), and determine its analytic structure. 
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The chromatic polynomial P{G, q) counts the number of ways that one can color a graph 
G with q colors such that no two adjacent vertices have the same color (for reviews, see 
I0I"[1])- The least positive integer q for which P{G,q) is nonzero is the chromatic number, 
x{G). Besides its intrinsic mathematical interest, the chromatic polynomial has an important 
connection with statistical mechanics since it is the zero-temperature partition function of 
the g-state Potts antiferromagnet (AF) [|5|, ^ on G: P{G, q) = Z{G, q,T = 0)paf- The Potts 
AF exhibits nonzero ground-state entropy 5*0 7^ (without frustration) for sufficiently large 
g on a given lattice graph and is thus an exception to the third law of thermodynamics. This 
is equivalent to a ground state degeneracy per site W > 1, since 5*0 = ksliaW. Denoting 
the number of vertices of G as n = v{G) and assuming a reasonable definition of {G} = 
lim„^oo G, we haveQ VFdG}, q) = lim„^oo PiG, qY^'^. Since P{G, q) is a polynomial, one can 
generalize q from Z4. to C The zeros of P{G, q) in the complex q plane are called chromatic 
zeros; a subset of these may form an accumulation set in the n — 00 limit, denoted B, 
which is the continuous locus of points where W{{G},q) is nonanalytic [|, |^ {B may 



be null, and W may also be nonanalytic at certain discrete points). The maximal region 
in the complex q plane to which one can analytically continue the function VFdG}, q) from 
physical values where there is nonzero ground state entropy is denoted The maximal 
value of q where B intersects the (positive) real axis is labelled qc{{G}). 

We consider strips of the square lattice with arbitrary length = m vertices and fixed 
width Ly vertices (with the longitudinal and transverse directions taken to be x and y). 
The chromatic polynomials for the cyclic and Mobius strip graphs of the square lattice were 
calculated for Ly = 2 in (see also 0|-[ll2|) and for Ly = 3 in |jl3|, After studies of the 



chromatic zeros for Ly = 2 in |P, [T^, |T6|, W and B were determined for this case in and 



for Ly = 3 in [T^. An important question concerns the effect of boundary conditions (BC's), 
and hence graph topology, on P, W, and B. We use the symbols FBCy and PBCy for free and 
periodic transverse boundary conditions and FBC^,, PBC3;, and TPBC^^^ for free, periodic, 
and twisted periodic longitudinal boundary conditions. The term "twisted" means that the 
longitudinal ends of the strip are identified with reversed orientation. These strip graphs 
can be embedded on surfaces with the following topologies:^ : (i) (FBCy,FBCs): strip; 
(ii) (PBCy,FBCx): cylindrical; (iii) (FBCj,,PBC2:): cylindrical (denoted cyclic here); (iv) 

^ At certain special points qs (typically qs = 0, 1, .., x(G')), one has the noncommutativity of limits 
liniq^q^ linin^oo ^'(G, q)^/" ^ hm„^oo liniq^g^ P{G,qY/^ , and hence it is necessary to specify the order of 
the limits in the definition of T4^({G}, qs) 0. We use the first order of limits here; this has the advantage of 
removing certain isolated discontinuities in W . 

^For the triangular lattice with cylindrical BC's, W and 13 were calculated in |0. Other calculations of 
P, W, and B have been performed for strips having BC's of type (i) jll-iol, (ii) (iii) [H-ll, (iv) 
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(FBCj^TPBa.): Mobius; (v) (PBCj^PSa): torus; and (vi) (PBCj^TPBa): Klein bottle.^ 
Here we present and analyze chromatic polynomials for the strip graph of the square lattice 
with Ly = 3 (i.e. cross sections forming triangles) and boundary conditions of type (v) and 
(vi): torus and Klein bottle. We recall that unlike graphs of type (i)-(v), the Klein bottle 
graph (vi) cannot be embedded without self-intersection in M^. For = m > 2 where they 
are well defined, the Ly = 3 torus and Klein bottle graphs have n = L^Ly vertices, e = 2n 
edges, the same girth g (length of minimum closed circuit) and number kg of circuits of 
length g, and the respective chromatic numbers X = 3 and x = 4. 

We label a particular type of strip graph as Gg and the specific graph of length = m 
repeated subgraph units, e.g. columns of squares in the case of the square strip, as {Gs)m- 
If one thinks of the graph as embedded on a rectangular strip of paper, with its upper and 
lower sides glued together and its longitudinal ends glued with direct or reversed orientation, 
then L3. is the length of this strip of paper in subgraph units. Writing 

n-l 

PmU,q) = Y.(-^yhn-jq''-' (1) 

j=0 



and using the results that ||25|, Q h„_j = (^^ ioi < j < g — 1 (whence /i„ = 1 and hn-i = e) 
and hn~{g-i) = {gli) ~ kg, it follows that for m greater than the above-mentioned minimal 
value, these /ij's are the same for the torus and Klein bottle of each type Gg. For a given 
Gs, as m increases, the hn-/s for the torus and Klein bottle graphs become equal for larger 
j- 

A generic form for chromatic polynomials for recursively defined families of graphs, of 
which strip graphs Gs are special cases, is 

P((a)™,?) = 5^c,(g)(A,(g))- (2) 

i=i 

where Cj{q) and the N\ terms Xj{q) depend on the type of strip graph Gs but are independent 
of m. 

For an Ly = 3, = m strip with (PBCj^FBC^.) one has P{sq{Ly = 3)^, PBCy, FBC^, q) 
q{q - l)(g - 2)(g^ - 6g2 + Uq - 13)""-\ whence 

W{sq{Ly = 3), PBGy, FBC, q) = {q' - 6q^ + Uq - 13)'/^ (3) 

with = 0. 



These BC's can all be implemented in a manner that is uniform in the length L^; the case (vii) 
(TPBCyjTPBCa;) with the topology of the projective plane requires different identifications as varies and 
will not be considered here. For connections between topology and graph properties, see e.g. pi, M. 
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In order to calculate P, one may use recursive methods based on iterative use of deletion- 
contraction theorems & |10|, or a coloring compatibility matrix method described in 



1^ , p6| . For the Ly = 3 torus (t) graphs, one finds 



Pisq{Ly = 3)^, PBCy, PBC, q) = Y1 ''tAK.r (4) 

i=i 



where 



Xt,i = -1 , ct,i = - + 8g - 1 , (5) 

At,2 = q'- Qq^ + Uq - 13 , q,2 = 1 , (6) 

(g-l)(g-2) 

At,3 = g - 1 , Q,3 = ^ , (7) 

V = g-4, Q,4 = (g-l)(g-2) , (8) 

At,5 = ? - 2 , Ci,5 = g(g - 3) , (9) 

At,6 = q-5, Q,6 = ^^^^ , (10) 

Xtj = -{q'-7q + 13) , ctj = q-l, (11) 

\8 = -(g-2)2, Q,8 = 2(g-1). (12) 
For the Ly = 3 Klein (K) bottle graphs one finds 

5 

P{sq{Ly = 3)„, PBCy, TPBC,, q) = J2 (13) 



where 



Xk,i = Xt,i = -1 , Ck,i = -{q - 1) , (14) 

Xk,2 = \t,2 = q^ - 6g' + Uq - 13 , Ck,2 = Q,2 = 1 , (15) 

X X 1 (g-l)(g-2) .... 

Ai^,3 = At,3 = q-1 , Ck,3 = -Ct,3 = , (16) 

AKA = At,6 = q-5 , ck,a = ctfi = , (17) 

Ax,5 = At,7 = -(g^ - 7g + 13) , ck,5 = Ctj = g - 1 . (18) 



The terms Xtj, j = 4, 5, 8 do not enter in eq. (|13]). We contrast this with earlier findings. For 
a given strip, Nx was found to be larger for (FBCj„PBC^) than (FBCj,,FBa) [||, |2l|, m. For 
the Ly = 3 square lattice strip case, Nx = 2 for (FBCy,FBC^) but Nx = I for (PBCy,FBC^), 
because of the special feature that the cross sections were complete graphs, Kp with p = 3 
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Figure 1: Chromatic zeros for the Ly — 3, 



10 torus graph. 



and hence the intersection theorem led to a factorized, monomial form for P. It was found 
||T3| , that for a given type of lattice strip, Nx is the same for the (FBCj/jPBC^;) = cyclic 
and (FBCy,TPBC^.) = Mobius topologies, although the Cj's were, in general, different. The 
present results show that reversal of orientation in the identification of opposite ends of a 
strip can lead to a change in A^a-Q 

Let C = J^f^i Cj- We find C = P{K^, q) = q{q - l)(g - 2) for the Ly = 3 torus graphs 



and C = for the L„ = 3 Klein bottle graphs. The zero results from the special constraints 



introduced by the boundary conditions and is analogous to the fact that C 
Ly = 2 Mobius square strip and its homeomorphic expansions 
Mobius strip graphs have C = 0; for example, for the Ly 
and kagome lattices, C = q{q 
Chromatic zeros for the L 



for the 
However, not all 



3 Mobius strips of the square 



q, respectively [P 



1) and C 

= 3, m = 10 torus graph are shown in Fig. |l]; with this value 
of m, the chromatic zeros for the Klein bottle graph are quite similar. The locus B and the W 
functions are the same for the torus and Klein bottle graph families. We find qc = 3, which, 
interestingly, is the same value as for the infinite 2D square lattice. The locus B has support 
for Re{q) > and separates the q plane into three regions. The outermost one, region Ri, 
extends to infinite |g| and includes the intervals q > 3 and g < on the real q axis. Region 



*A different sort of change in P, accompanied by a change in B, can be obtained if one considers a 
homogeneous recursive family and the same family with a finite inhomogeneous subgraph inserted, e.g., the 



■'rope ladder" graphs of ]lq | or two such subgraphs forming ends, viz., the J{Y[H)I strip graphs in |18 
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i?2 includes the real interval 2 < q < 3 and extends upward and downward to the complex 
conjugate triple points on B at qt and q^, where qt — 2.5 + lAi. Region is the innermost 
one and includes the real interval < q < 2. The boundary between R2 and R3 curves to the 
right as one increases \Im{q)\, extending from q = 2 upward to qt and downward to q^. As 
is evident in Fig. |l], the density of chromatic zeros on the Ri — R3 boundary near g = and 
on the i?3 — i?2 boundary is somewhat smaller than on the right-hand part of the R2 — -Ri 
boundary. 

In region Ri, Xt^2 = ^k,2 is the dominant Xj, so 

W ={q^- 6q^ + Uq - 13)^/^ , q e Ri . (19) 

The fact that this is the same as W for the (PBCy,FBCj,.) case, eq. (H), is a general result. 
The importance of the PBCy is evident from the fact that for the same width of three squares, 



the strip with (FBCy,FBC^) yields a different W p 
In region R2 Xtfi = Xk,4, is dominant, so 

\W\ = \q-5\^/\ qeR2 (20) 

(in regions other than only \W\ can be determined unambiguously [0). In region i?3, 
= Ax,5 is dominant, so 

\W\ = \q^ -7q + l3\^'\ qeR3. (21) 

The outer boundary separating Ri from the inner two regions is oblate, extending out to a 
maximum of about \Im{q)\ ^ 2.5 for Re{q) ~ 1.5 (and passing through q = and 3). All 
the three points, q = 0,2,3, where B crosses the real q axis, it does so vertically. The present 



results are in accord with the inference [13, O] that for a recursive graph with regular lattice 



structure, a necessary and sufficient condition for B to separate the q plane into two or more 
regions is that it contains a global circuit, i.e. a path along a lattice direction whose length 
goes to infinity a.s n oo; here this is equivalent to PBC^;. The fact that B is the same for 
these torus and Klein families means that none of Xtj, j = 4, 5, 8 is a dominant term. 

Our calculations of the zero-temperature Potts antiferromagnet partition functions (chro- 
matic polynomials) and exponential of the entropy, W, for Ly = 3 strips of the square lattice 
with periodic transverse and periodic and twisted periodic longitudinal boundary conditions 
(torus and Klein bottle graphs) thus elucidates the role that these boundary conditions and 
the associated topologies play; the torus and Klein bottle graphs have interestingly different 
chromatic polynomials, with different A^a, but the W functions and hence the boundaries B 
are the same. 

The research of R. S. was supported in part by the U. S. NSF grant PHY-97-22101. 
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